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Secure Transmission in Amplify-and-Forward 
Diamond Networks with a Single Eavesdropper 

Siddhartha Sarma*, Samar Agnihotri^ and Joy Kuri* 


Abstract 

Unicast communication over a network of M-parallel relays in the presence of an eavesdropper 
is considered. The relay nodes, operating under individual power constraints, amplify and forward the 
signals received at their inputs. In this scenario, the problem of the maximum secrecy rate achievable 
with AF relaying is addressed. Previous work on this problem provides iterative algorithms based 
on semidefinite relaxation. However, those algorithms result in suboptimal performance without any 
performance and convergence guarantees. We address this problem for three specific network models, 
with real-valued channel gains. We propose a novel transformation that leads to convex optimization 
problems. Our analysis leads to (i) a polynomial-time algorithm to compute the optimal secure AF rate 
for two of the models and (ii) a closed-form expression for the optimal secure rate for the other. 


1. Introduction 

Wireless communication, by its inherent broadcast nature, is vulnerable to eavesdropping by 
illegitimate receivers within communication range of the source. Wyner in [[^, for the first 
time, information-theoretically addressed the problem of secure communication in the presence 
of an eavesdropper and showed that secure communication is possible if the eavesdropper 
channel is a degraded version of the destination channel. The rate at which information can be 
transferred from the source to the intended destination while ensuring complete equivocation 
at the eavesdropper is termed as secrecy rate and its maximum over all input probability 
distributions is defined as the secrecy capacity of channel. Later, Q extended Wyner’s result to 
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Gaussian channels. These results are further extended to various models such as multi-antenna 
systems 0. 0 . multiuser seenarios 0, fading ehannels Q, 0- 

An interesting direction of work on secure eommunieation in the presence of eavesdropper(s) 
is one in which the source communieates with the destination via relay nodes [[9|-pA|. Sueh work 
has eonsidered various seenarios sueh as different relaying sehemes (amplify-and-forward and 
decode-and-forward), constraints on total or individual relay power eonsumption, one or more 
eavesdroppers. However, tight eharaeterization of seerecy capaeity or even optimal aehievable 
rate is not available except for a few specific scenarios. 

We eonsider unicast communication over a network of M parallel relays in the presence of an 
eavesdropper. The relay nodes, operating under individual power eonstraints, amplify and forward 
signals reeeived at their inputs. The objeetive is to maximize the rate of seeure transmission from 
the source to the destination by choosing the optimal set of scaling factors for the AF-relays. 
In general, the problem is non-eonvex and global optima cannot be guaranteed. However, we 
establish that convexification of this optimization problem is possible for degraded eavesdropper 
channels scenario where each relay to destination ehannel gain is larger than the eorresponding 
relay to eavesdropper ehannel gain. Though sueh a scenario is limited, we argue that it may 
still provide insights into the nature of the optimal solution for general problem and help us 
construct low-complexity eapacity achieving achievability sehemes. Later, we diseuss two special 
scenarios: scaled eavesdropper channels and symmetric network and provide a polynomial-time 
algorithm and a elosed-form solution, respectively, to eompute the optimal seeure AF rate. 

To the best of our knowledge, we are the first to propose a eonvex optimization based solution 
to obtain the global optimum for the degraded ehannel seenario, where channel gains are real 
number^ Previously, the problem addressed in this paper was also eonsidered for general ehannel 
seenarios by pT| for a single eavesdropper and pA| for multiple eavesdroppers. However, the 
semidefinite relaxation based approaches proposed therein only provide suboptimal performanee 
without any guarantee on the performance (even for simple network models) and the convergence 
of the proposed iterative algorithms. Unlike sueh previous work, we propose a novel variable 
transformation that results in eonvex optimization for the degraded ehannel seenario. 

Our main eontributions are as follows: 

'The more general case with complex channel gains is currently under investigation. 
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Fig. 1. The diamond M-relay AF network with an eavesdropper 


• For the degraded eavesdropper ehannel model, we propose a novel transformation for the 
non-convex problem and present an iterative algorithm with guaranteed eonvergenee to 
eompute the optimal seereey rate aehievable with AF relaying. We also provide a lower- 
bound to the optimal rate based on the zero-forcing approaeh and an upper bound based on 
the total relay power eonstraint. 

• For the network model with sealed eavesdropper ehannels, motivated by the greedy ap- 


proaehes proposed in [14|- [15| to eompute the optimal AF rate over diamond networks, 
we propose a novel efficient scheme to compute the optimal secrecy rate achievable with 
AF relaying. 

• Finally, for the symmetric network model, we provide a closed form expression for the 
optimal secrecy rate achievable with AF relaying, along with the corresponding optimal 
vector of scaling factors. 

This paper provides some of the key results of our work. Various extensions and detailed 


discussions can be found in our ArXiv submission p^ . 

Organization: Section |I^ introduces the system model and notation. In Section we analyze 
the optimal AF secrecy rate for three specific network models. The performance of proposed 
approaches is evaluated in Section Finally, Section |V] concludes the paper. 


II. System Model 

Consider the Gaussian “diamond” or parallel M-relay network as depicted in Figure The 
source node s transmits a message to the destination t with the help of M parallel relays. 
However, the signals transmitted by the relays for the destination are also overheard by the 
eavesdropper e. The channel inputs at time n e N at the source s and the relay are denoted 
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by Xs\n\ and Xj[n], respectively. The channel outputs at time n e N at the destination t, the 
eavesdropper e, and the relay are denoted by yt[n],ye[n], and yi[n]. These channel inputs 
and outputs are related as: 


yi[n] = hsiXs[n] + Zi[n] 


M 


yk[n] = Yj hikXi[n] + Zk[n], ke {t, e} 


( 1 ) 

( 2 ) 


2=1 


where {zi[n])i^n, izt[n])n, and {ze[n])n are independent and identically distributed Gaussian 
random variables with zero mean and variance independent of channel inputs. The channel 


gains {hsi)fLi, and {hie)fLi are real numbers, constant over time ( [ 151, [ 17|) and known 

(even for the eavesdropper channels) throughout the network pO|-p^. 

In amplify-and-forward relaying, each relay scales its received signal before transmitting. The 
maximum scaling factor at a relay is determined by its individual power constraint. We consider 
the following average transmit power constraint: 

£'[a;^[n]] < P,, —oo < n < oo, z e (s, 1, 2, ..., M] 

Therefore, the channel input at the z*^ relay can be written as: 


Xi[n + 1] = Al/iN, 


(3) 


where the scaling factor (3 is subject to the constraint 

0 < (4) 

Assuming equal delay along each path, for the network in Figure the copies of the source 
signal (Xs[.]) and noise signals (^i[.]), respectively, arrive at the destination and the eavesdropper 
along multiple paths of the same delay. Therefore, the signals received at the destination and 
eavesdropper are free from intersymbol interference (ISI). Thus, we can omit the time indices 
and use equations ([T])-@ to write: 

M M 

Vk hgif3ihi}iXs + flihij^Zi + 2 :^, k g {f, e} (5) 

2=1 2=1 

The secrecy rate at the destination for such a network model can be written as 0, R.(P,) = 
[I{xs; yd) — I{xs] where I{xs] y) represents the mutual information between random vari¬ 

able Xs and y and [zz]+ = max{zz,0}. The secrecy capacity is attained for Gaussian channels 
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with the Gaussian input ~ J\f{0,Ps), where E[x^] = Pg, Q. Therefore, the optimal secreey 
rate aehievable with AF relaying ean be written as the following optimization problem. 


Rs{Ps) = max [Rd{Ps, f3) - Re{Ps, /3)] 
ri l + SNRdiPs,(3) 


where SNRk{Ps,f3) = 'ygTk{f3), k e {t,e} with 


r.(/ 3 ) = 


(hs,k‘/3)^ 


l+/3*Dk/3’ 

7s = K,k = [hsihk, hs 2 h 2 k, • • • , hgMhukW and Dk = diag{[hlf,, 

diagonal matrix with diagonal entries written in veetor form. 

In the following section, we analyze the optimal secure AF rate problem in 
specific network models. 


(6a) 

(6b) 

(7) 

,hMk]) isaMxM 


for three 


III. The Optimal AF Secrecy Rate Analysis 

We discuss the optimal solution of the maximum AF secrecy rate problem @ in the M-relay 
diamond networks with single eavesdropper in degraded eavesdropper channels scenario. Later 
we consider two special cases, namely scaled eavesdropper channels and symmetric network. 


A. Degraded Eavesdropper Channels 

Consider a general channel model where gains along the source to relays and relays to the 
destination channels may take arbitrary values. However, the relay to eavesdropper channels are 
degraded versions of the relay to destination channels, that is, hu > hie or he = cx^ht, ai < 
l,ze {!,...,M}. 


1) Optimal secrecy rate for individual relay constraints: From equation @ we can rewrite 
an equivalent optimization problem in the following manner: 

1 + SNRM 


max max 
V /3 


I + r] 

subject to: SNRe{f3) ^ rj, 


(8a) 

(8b) 




(8c) 
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The optimization problem @ can be solved in two steps. At first, we fix an rj and solve the inner 
optimization problem for optimal /3. We first present an approach to solve the inner optimization 
problem, then we address the second step involving the variable rj. 

We perform the following transformation on the inner optimization @ to convert the non- 


convex problem into a convex one. Vi = Ui/{1 + ob equivalently u) = 

The objective function becomes v^hghgV and it can be shown that it suffices to maximize the 
linear objective h^v. The reformulated inner optimization problem is: 


Vl—"v 


max hgV 

V 

subject to: v*Ce(r 7 )v ^ 1, 

v*DiV ^1, i e {1, 2,... M}, 


(9a) 

(9b) 

(9c) 


where Ce{v) = 


heh* Pg 
r] 0-2 


+ I-De, he = [hsipi,hs2p2, ■■■,hsMpMf, De = diag{[pl, p^,..., pj^]) 


and pi = ^, VT Di characterizes the constraint for each relay node and is given by 


1 + 




if m = n = i 




^id^i^max 

1, if m = n A * 


0, otherwise 

For the degraded channels pi < 1, so Ce{p) is a positive semidefinite matrix. As the objective 
function is linear and constraints are convex, so the problem in @ is convex and can be solved 
optimally using well-known numerical techniques e.g. interior-point methods. 

Once we have obtained the solution for the inner optimization problem over transformed 
variable v, then ( [T^ reduces to a one dimensional optimization problem. If we denote the 
optimal solution for the inner optimization problem by v*(? 7 ) for a fixed p, then the following 
proposition describes the nature of the objection function ( [Sa] ). 

Proposition 1: The objective function f{p) = is unimodal in p for p e [0, oo). 

Proof: Please refer to Appendix ■ 

As fir]) is unimodal, so we can guarantee the convergence to global optima using line search 
algorithms {e.g. golden section [18|). Instead of searching over whole positive real line (M+), 
we can limit the range of p by further analysis. 
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Range of r]: As SNR is always positive, so r] = 0 serves as a lower bound. This rj value 


eorresponds to the zero-foreing solution (III-A31 whieh serves as a lower bound for optimization 
problem ([^. For upper bound on rj, we solve the following optimization problem: 

r]max = maxTeif3), subjeet to: /3*A/3 ^ Ptot = Pi+ ■■■ + Pm 

where A = diag[hl;^Ps + • • • , h^j^Pg -1- cr^]. The eonstraint considered here is total 

power constraint and corresponds to a larger feasible solution space as compared to individual 
constraints. Therefore, T]max calculated here is an over estimate of the maximum value of t] 
with individual constraints. The solution of this problem is rjmax = hs e*De where De = 

Along with the solution of the optimal AF secrecy rate in this setting, we also provide 
upper and lower bounds on the optimal rate. For upper bound, we replace the individual power 
constraint on each relay by a single total power constraint on /3 vector. For lower bound, we 
consider the zero forcing solution d). [ [T3l , where /3 values are chosen such that the transmitted 
signals get canceled at eavesdropper. 

2 ) Optimal secrecy rate for the total power constraint: From equation ([^ the total power 
constraint optimization problem can be expressed in following manner: 

1 + SNRM 


max max 

7? /3 


1 + p 

subject to: SNRe{f3) < p, 
/3*A/3 < Ptoti 


(10a) 

(10b) 

(10c) 


One can check that the total power constraint ( |10c| ) results in larger feasible solution space 
as compared to (|^. Therefore, the optimal objective function value obtained will upper bound 
the original individual constraint problem ([^. Using the transformation described above, we can 
reformulate the optimization problem in the following manner: 


max hgV (Ha) 

V 

subject to: v*Ce(r 7 )v ^ 1, (Hb) 

v*DtV ^ 1 


(11c) 








where Dt = diag{[(iX 2 , ■' ’ )Cm]), 0 = (1 + )• As Dt is a diagonal matrix with 

positive diagonal elements, so ( |lle| ) is a eonvex eonstraint. This problem ean be solved in a 
manner similar to the one for the individual power eonstraint. 

3) Zero forcing for individual constraints (rj = 0): In this approaeh we equate SNRe to 
zero while maximizing SNRd. By imposing this extra eonstraint we are redueing the feasible 
solution spaee, therefore, the solution obtained will provide a lower bound to original problem. 
The optimal seeure AF rate problem for this ease ean be written as: 

P [Tiii hsiPihit 


max 

/3 




M 


subjeet to ^ hsi^ihi^ = 0 , 


i=l 


/di^max ^ fdi ^ 


( 12 a) 

( 12 b) 

( 12 e) 


We reformulate it as the following quadratie program. 

Proposition 2: The problem in ( [T^ is equivalent to the following quadratie program whieh 
ean be solved effieiently 

max w*w, subjeet to: Hw = [1 0]\ H/ 3 W < 0 

W 

Proof: Please refer to Appendix ■ 

Remark 1: As this formulation does not rely on the degradedness assumption, it ean be used for 
general ehannels also. 


B. Scaled Eavesdropper Channels 

Sealed eavesdropper seenario is a speeial ease of degraded ehannel seenario. Here, all the 
elements of veetor ol have same value, i.e. ai = a < l,i e Therefore, eavesdropper 

ehannel ean be expressed as = aht. Though the approaeh diseussed in the last seetion ean 
be used to solve the optimal AF seereey rate problem eorresponding to this network model, 
here we propose another seheme that better exploits the eorresponding ehannel properties, thus 
resulting in a more effieient implementation. Before going to details, we provide an overview 
of the proposed solution to the sealed eavesdropper problem. 






9 


1) Due to the structure of the optimization problem we identify that a common function is 
present in both numerator and denominator and maximizing that function will result in 
maximization of the original problem. 

2) Consider a new variable uji = hidPi, its upper bound Ui^max = hidPi^max, and a new 

parameter gsi = e {1,2, Now we relaxed the derived optimization 

problem by replacing all the individual constraint by a single total constraint which can 
be solved efficiently. 


M 


Original: max f{uj) = ^ gs^iUi subject to: Ui ^ uji^„^ax i e (1,2,..., M} 


' M 


M 


Relaxed: max j ^ gs,i<^i subject to: ^ 


2 2 
uf = r 




2=1 


Here, e [0, ||<-c;|p] 

3) If the solution of this relaxed problem satisfy individual constraints, then we have solved 
the problem optimally. 

4) If the solution of that total constraint problem violates individual constraints, then we have 
to impose the individual constraint criteria to obtain the solution. In other words, some 
of the variables have reached their boundary value and others are still inside constraint 
region. The question is how to identify the variables that have reached their respective 
upper bounds? 

Let us begin with the obvious scenario, where = lluJmaxlP- The solution is trivial i.e. 

V{ 1,2,..., M}. Now, if wc Want to reduce <5 from all the variables, how 
should we proceed so that objective is least effected. Consider, 0 < < 1 and ^ = 1. 

i 

If we reduce 59i from each variable, then for 5-^0, neglecting 5^, the rate of change in 
objective function with respect to can be written as: 

9s,i 


^ = y = V/9. 

(ir2 ^ dui dojj ^ 




u. 


9s,i* , . 9s,i 

-where i = arg mm 


i,max 


,max 


00 . 


i,max 


Therefore, we should reduce the variable cuj*, to keep the rate of decrement of the ob- 
jcctivc function minimum. AA/c c3-n keep reducing only till -y- = —, where 


CJ'5}C 


9s,j 


J" = arg . From that point onwards variables ooi* and ooj* both has least 

effect on objective function among all variables and, therefore, we keep reducing them 
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jointly till we reach 


9s,k^ 


where k* = argmin^^ 


9s,k 
UJk rna: 


and it continues in similar 


,max 

manner. 

From the above discussion one can conclude that an ordering among the variables are pos¬ 
sible depending on their contribution towards objective function. Variables that contributes 
more towards objective function have higher tendency of violating individual constraints 
as compared to other variables. So, we sort the variables in descending order based on 
there value and re-index them in following manner: 

^i,Tnax 

9s, {1) ^ 9s, (2) ^ ^ 9s,{M) 








^(1) ^max 01(2) ,max ^{M),max 

Due to re-indexing, we can say that a;(^) has higher tendency of violating individual 
constraint than Now, we can divide the range [0, ||<-Umax|p] based on this ordering. 

An interval [r^, is defined in following manner: 

,2 M 


= 

' m 


i=l 




^{i),max 


+ 


(m) ,max 


9 


Tj 3. 


s,{m) i=m+l 


2 

's,{i) 


In simple words, within this interval variables indexed by (1) to (m) has reached their 
upper bounds and we need to solve the total constraint problem for rest of the variables. 

5) We argue that the solution for each of these subproblems within an interval can be obtained 
by solving a one-dimensional maximization problem. 

6) We choose the solution that maximizes the original objective function over those sub¬ 
intervals. 

We rewrite the achievable AF secrecy rate in terms of cu: 


Rs{u;) = ^log 


1 -I- 


1 -f 

where = (j]fii9si<^i^ , Qi = and g 2 = with So, in this setting 

the problem in @ is: 

Rs{Ps) = max Rs{oj) (13) 


In the rest of this subsection we provide an approach to efficiently compute the optimal 
solution of the problem in ( [T3] ) and then summarize it in an algorithm. 

Note that a non-zero secrecy rate can be obtained only if a < 1. Thus, > ^ 2 - This implies 
that is maximized when T'(lj) is correspondingly maximized. To maximize T'(lj), we 















11 


formulate the following problem: 


M 


M 


2 2 
oji = r 


(14) 


max I ^ (y'sjcej j subjeet to ^ 

V*=l / i=l 

It is easy to prove that the optimal a;(r) veetor for this problem is equal to u;(r) = By 

replaeing this in the seereey rate equation ( [T3l ) along with the sum eonstraint we formulate the 
following optimization problem in r: 


1 + 


max 

r 




1 + 


(15) 


1+a^r 


Lemma 1: The optimal solution to the problem in ( [T5] ) is r* = a^Jl + ||gs|P)^^ 

Proof: By differentiating the objeetive funetion of ( [T5| ) with respeet to r and then it equating 
to 0, we get the solution. ■ 

Now, if satisfy individual eonstraints of ce, then this is the optimal solution, else we 

must eonsider the next step that eonsiders seenarios where the solution obtained using above 
method violates any of the aij’s eonstraints. 

Arrange aij’s aeeording to and denote them as a;(i) ^ ^ u]{m)- The gsi are ordered 


aeeordingly with the ordered values denoted as gs,{i)- Then, we divide the feasible set of the 
problem in ( [T?] ) in M non-overlapping subsets, [r^, r^+i], m e {0,..., M — 1}. The interval that 
results in the largest value of the problem in ( fTS] ) solved over it as diseussed below, is ehosen 
as the one that results in the optimal solution of ( [T3] ) over all intervals and the eorresponding 
optimal veetor of sealing faetors /3*. The eorreetness of this approaeh follows from an argument 
similar to the one in [ [T4l Lemma 1]. 


For the ease of exposition, define Pm = 9s,{i)^{i), 

yM 2 

Z-J2=m+1 


max^ Hm 


a = 

qm Zj2=1 




, 3.11(1 


The intervals can be calculated as follows. The first interval is [0, ri], where ri corresponds to 
the value at which the first variable (after ordering) is equal to its upper bound. Therefore, cj* = 

= jg|ri. Thus, ri = = yyyyryh- 

Similarly, = ^{sm/ me{l,... ,M- 1} and tm = qu- 

Remark 2.- If 0 ^ r ^ ri, then we ean obtain the optimal sealing veetor using the solutions of 


((U]) and ([T5]). 

Consider evaluating ( [T3| ) over the interval [r^, rm+i]. In this stage the first m ordered variables 
have reaehed their respeetive boundaries. These variables are replaeed by their respeetive upper 
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bounds. For the rest of the variables we search for their respective optimum values as discussed 
next. 

The part of the objective function in ( fT?) ) depending on variables \uj(^rn+i), ■ ■ ■ ,i^(m)] can be 
upper bounded as: 


^ Am 11gs,{(m+l),-,(M)} | 
for all ||u;{(m+i),.,(M)}|| = - Qm, 


where U{(m+i),.,M} := [u(m+i), ■ ■ • , W(m)]* and Am > 0. In other words, the objective function is 
maximized when the vector variables [ci;(m+i), • • • > lies in the direction of [gs,{m+i) : 9s,( m)]^ 

vector. Therefore, we have to find the optimal scaling factor Am for those variables. Replacing 
the first m ordered variables by their respective upper bounds and rest of them by a scaled vector 
of gs,{(m+i),-,(M)}, we obtain the following optimization problem in terms of Am- 


max 


1 I ^ (.Pm ~t~ ^ 

a^{Pm + Sm>^m)‘^ 


\m-.\m>0 1 + ^ a (Pm + SmAm) 

Is 1+Q,2q^+Q,25^_\2 

The solution to the above problem can be calculated by solving the following degree-4 
polynomial equation: 


(1 gm)(l + a'^qm) _ Pm(7sPm«^ + 2gm«^ + a^ + 1) 

S^a2(l 7^Sm) + 'IsSm) 

_2 Pm(2 -+■ 37gSm) .^3 _ \4 _ ri 

The coefficients of this equation are all negative except the first one. Using the Descartes’ rule 
of sign there is only one positive root as there is only one variation in sign of the coefficients. 
Given that 0 < Am, it is the desired solution. Note that, if Am > W{m+i),max/gs,{m+i), then the 
solution will not belong to predefined sub-interval (•.• = (sm+i/5'J(m+i))'^(m+i),max+'?m+i < 


Am-Sm+rf Qm+l)- 


We summarize the whole procedure in Algorithm 


C. Symmetric Network 

A special case of the general s — t and s — e diamond networks introduced above is the 
symmetric setting, where hsi = hs,hit = ht, hie = he and Pi = PR,i e M}. In this 
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Algorithm 1 Algorithm for Scaled Eavesdropper Channel 

Input* O', Qs, ^i,maxi Pi ^ 

Output: uj* 

1 : Calculate r* = , / — and evaluate uj = 

yaVl+l|gs||2 

2: if (jJi ^ ^i,maxi ^ ^ {I7 • • • 1 thCU 

3 : (jj* = (jj (Lemma [T]) and exit. 

4: else 

5: Sort cCiS according to —2^^ in the descending order. 

^i,max 

6: for m = 1 to M — 1 do 

7: Solve polynomial equation ( [T^ to calculate the A^- 

8: Evaluate cUm using Am corresponding to interval [rm,rm+i] starting from first i for 

which Xi < W(^i+i)^rnax/gs,ii+i)- Dcnotc that i as zq. 

9 : end for 

10 : (jj* = (JL>^, (jj^ maximizes the solution of ( fT3| ) over all cUm, m e (zo,..., M — 1}. 

11 : end if 


setting, for a given vector of scaling factors, /3 = {(3i,..., (3 m), the achievable AE secrecy rate 
is: 


Rsi(3) = ^log 


( 


1 + 


M 


\ 


\ 


7I S A 

\ii=i 
M 

^+2 A ? 
1=1 / 


( 


1+ 


V 


!i(M^ 

M 

PI 

1=1 / J 


where % = 'fghs, u = and g, = Note that Rs{(3) is non-zero only if z/ < /z, i.e. ht > hg. 
Thus, using the optimal secure AE rate problem in this setting is: 


Rs{Ps) = Rs{(3) (17) 

Proposition 3: Eor the symmetric M-relay diamond network, the optimum scaling factors for 
all relays are equal, i.e. /?* = = • • • = = (3*. 

Proof: Please refer to Appendix ■ 

Corollary 1: Eor symmetric M relay diamond network, the optimal value of the scaling factor 
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for each relay is: 


= U 


a 


PuM^hm 


1/4 


t^max 5 


(J 


Pr(3^. 


2 

max 




f^maxi O.W. 

Proof: Using Proposition the optimal AF secrecy rate problem in @ is reduced in this 
case to the following single parameter optimization problem: 

Ps \ 


max 


1 + 


P, 


1 + 


a2 1 + M/32/^2 j j ■ ^2 1 + M/ 52 /i 2 j 

With the parameters: = -s/Mhg, Kd = -sfMht, and hre = VMhg, this problem reduces to 

the single relay AF secrecy rate maximization problem, addressed in Lemma below. ■ 

Lemma 2: The optimal AF scaling factor for single AF relay network is: 


(^AF = •* 


a 


PRhLhh 


1/4 


t^max") 


a 




2 

max 


< hzA, 


rd ^re 


R ^rd ^re _ 

Pmax") O.W. 

Proof: Please refer to Appendix ■ 

Remark 3: One can obtain an upper bound on secrecy rate for degraded channel scenario 
( III-A[ ) by using the result of symmetric channel scenario. For example, if we form a M relay 
symmetric network by considering hg = maxjhs}, ht = max{ht} and he = min{he}, where 
max{x} and min{x} denotes the maximum and minimum element of vector x, respectively. 


IV. Results 


To evaluate the performance of the scheme proposed in subsection III-A consider a network 
whose main channel gains are sampled from a Rayleigh distribution with parameter 0.5. The 
degraded channels for eavesdroppers is obtained by multiplying a factor (~ Uniform[0, 1]) with 
the relay to destination channel gain. We average the results of 1000 such network instances 
while plotting the graphs. 

In Figure we plot the variation of the optimal AF secrecy rate with respect to the source 
power (Pg) for the degraded eavesdropper channels. The zero-forcing lower bound and the total 
power constraint upper bounds are also plotted. 
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Source Power (P^) dB 


Fig. 2. Plot of the optimal secrecy rate (Rf) versus source power (Ps) for five relay node diamond network using the scheme 

2 


proposed in subsection 


III-A 


Flere we consider Pr = 5 and cr = 1. 


V. Conclusion 

The problem of maximizing the rate of seeure unicast transmission over a network of parallel 
AF relays in the presence of an eavesdropper is considered. Unlike previous work that provides 
iterative algorithms with suboptimal performance without any performance and convergence 
guarantees, our work provides either a closed-form expression or polynomial-time algorithm to 
compute the optimal secure AF rate for three specific network models considered. In future, we 
plan to investigate optimum secrecy rate of general AF networks. 
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Appendix A 

Proof of Proposition [I] 


Proof: In optimization problem (|9]), the individual eonstraints ( |9c| ) are independent of rj. If a 
vector xi satisfies the constraint ( |9bl ) for rji, then xi will satisfy it for t ]2 > rji. 1 x^Cg(r)i)xi = 
^ + ^iDeXi > + x^DeXi = x*Cg(r) 2 )xi, here De = I - D^. In fact, as i] 

increases feasible set for ( |9b| ) also increases. As a result the feasible set for optimization problem 
(|^ also increases. Therefore, (hgV*(7)i))^ ^ (hgV*(7)2))^- Recall that, t] = 0 corresponds to zero¬ 


forcing solution (Sec. Ill-A3) and, therefore, /(r)) increase with rj initially. But later constraint 


( |1 lc| ) become active and constrain the feasible set. As a result, the numerator of f{r]) remains 
fixed, whereas, denominator increases with r]. This results in decrease of /{t]). ■ 
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Appendix B 

Proof of Proposition [2] 

Proof: It can be shown that the optimal solution does not change if we rewrite the objective 


function as min because id A 0. If we consider a new variable v 

(2i=l hsi^ihid) 


such that hsi(3ihid = v, then in terms of variable vector w = [ 7 *, we can easily 


write a quadratic optimization problem. The denominator of the objective function ( |12a[ ) can 
be written as a constraint in terms of this variable vector as, [h*,0]w = 1. The eavesdropper 
constraint can be expressed as [h*, 0]w = 0. We combine them as matrix H and write as equality 
constraint. H/j can be obtained by rearranging the following constraint: —hd,il3i,maxWM+i < w* ^ 


Appendix C 

Proof of Proposition [3] 
Proof: Taking derivative of Rg with respect to f5i, we get 


M M 

(2 _ fl.Y V 




i=l 

M 


i=l 


M 


M 
>7 = 1 


i=i 

i=l 


M 




0 Eft 

1 

II 

0 

7=1 

J 7=1 


optimal solution and fr ¥= 0, 

hence 


^ M / M \\ / M \ / M 

Sft PS/’A-' + a' - E/S? 


y,2 = l 


0=1 


i=l 


\ 2=1 


M 

( Zj a ) + 1^/1 = 0 (18) 

0=1 


Similar equation can be obtained if we take derivative with respect to /?,, j ¥= i and j e 


{1,2,.. .,M}. Subtracting this equation from equation (18 1 we get 


M 


M 


M 




= 0 


0 = 1 


V7 = l 


Hence, /?* = /3j, i,j e {1, 2,..., M} 


(19) 





Appendix D 
Proof of Lemma [2] 
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For single relay network the problem is: 

^ / P, hll3^hl^ 

13 af = argmax 1 + 

/3 V 


1 + 


P. hlP^hl, 

(T^ 1 + 


1 + 

Differentiating the objective function with respect to (3 and setting the derivative equal to zero, 
we get for the optimum f3 value 




a 


A 


2 

max 


LPRh'X,^ 

Further, given that the maximum value that (3 can take is bounded from above by (3max and the 
mono tonic increasing nature of the objective function from /3 = 0 to /3 = {3 max-, we obtain the 
desired result. 







